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I.  INTRODUCTION 


There  was  a  need  for  the  equation  of  state  of  porous  propellant  so  that 
code  calculations  could  be  made  to  predict  the  response  of  propellant  beds  to 
shock  pressures  of  varying  intensity  and  duration. 

The  work  reported  here  is  related  to  the  XM1  Compartmentalization  Study. 
One  threat  to  crew  and  tank  is  propellant  initiation  caused  by  shaped  charge 
jet  penetration.  What  shock  pressure  conditions  will  be  produced  in  the 
propellant  bed  during  shaped  charge  penetration?  Will  the  propellant  reaction 
under  these  conditions  be  so  violent  that  quick  venting  is  not  possible?  In 
order  to  answer  these  questions,  it  is  helpful  to  know  the  shock  Hugoniot  of 
the  impacted  material.  The  work  reported  here  was  specifically  directed 
toward  the  measurement  of  the  shock  Hugoniots  for  live  and  inert  propellants 
and  calculation  of  the  shock  Hugoniots  of  porous  beds  of  these  materials.  The 
live  propellant  material  is  the  seven  perforation  M-30  propellant  which  is 
used  in  the  105-mm  SPDS  and  HEAT  rounds.  The  inert  propellant  consists  of  a 
thermoplastic  material  which  has  been  extruded  to  the  same  configuration  as 
the  live  propellant  and  has  approximately  the  same  density^  The  porous 
Hugoniots  were  calculated  from  the  solid  Hugoniots  U3ing'"the  Mie-Gruneisen 
equation. 


II.  SAMPLE  PREPARATION 


A.  M-30  Propellant 

The  propellant  come3  in  the  form  of  a  grain  approximately  1.59  cm  long, 
0.71  cm  diameter  with  seven  cylindrical  perforations  0.08  cm  diameter  within 
the  grain,  parallel -to  the  long  axis,  as  shown  in  Figure  1.  Samples 
approximately  0.32  cm  long  were  machined  from  several  propellant  grains  and 
carefully  sanded  and  measured  for  thickness.  The  small  sample  thickness  wa3 
necessitated  by  the  small  diameter  of  the  propellant  and  the  need  to  make  a 
shock  velocity  measurement  through  the  sample  before  rarefactions  originating 
at  the  sample  diameter  boundary  can  influence  the  shock  wave.  The 
perforations  were  filled  in  with  modeling  clay  in  order  to  prevent  jetting  and 
eliminate  rarefaction  effects  on  the  shock  wave  in  the  sample.  The  clay  had  a 

3  3 

density  of  1.65  gm/cm  ;  the  propellant  density  was  1.66  gm/crn  ,  M-30  is  a 

tripj.v.  iase  propellant  which  has  the  .  follow! ng  ingredients: 


nitrocellulose 

28.0 

f 

nitroglycerine 

22.5 

% 

nitroguanadine 

47.7 

% 

ethyl  centralite 

1.5 

% 

cryolite 

0.3 

% 

B .  .Inert  Simulator 

The  inert  simulator  is  an  extruded  thermoplastic  made  by  Radford  Army 
Ammunition  Plant  and  it  was  made  to  simulate  the  density  of  M-30  propellant. 
It  consists  of  the  following: 
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7  CYLINDRICAL  HOLES 
HAVING  A  DIAMETER 
OF  .08  cm 


Figure  1.  Initial  Configuration  of  the  M-30  Propellant  Grain. 


cellulose  acetate  64 .2  % 
triacetin  (plasticizer)  15.0  % 
graphite  2.4  $ 
red  lead  ( Pb^O^ )  1  18.4  % 


The  inert  simulator  has  a  density  of  1.55  gm/cm  .  It  comes  in  the  form  o 
pellets  having  approximately  the  same  dimensions  and  configuration  as  the 
live  propellant.  In  order  to  obtain  uniform  samples  suitable  for  our 

measurements  the  inert  was  heated  to  100°C  and  compressed  in  a  metallurgical 
press  at  10,000  psi  to  get  a  solid  cylinder  3. ,18  cm  diameter  by  3.81  cm  high. 
Samples  .95  cm  thick  were  machined  from  this  and  measured  for  thickness. 


III.  PROCEDURE 

The  impedance  matching  technique1  was  used  to  determine  the  solid 
Hugoniots  for  the  live  and  inert  propellant,  materials.  We  used  a  -l0*'16  e“ 
diameter  plane  wave  lens,  a  2.54  cm  thick  TNT  pad,  and  various  buffer  plates 
to  provide  a  range  of  known  pressures  in  the  buffer  plate.  The  sample  of 
known  thickness  was  placed  on  the  buffer  plate  and  front  lighted  by  an  argon 
bomb  light  source.  The  transit  time  of  the  shock  through  the  sample  was 
recorded  by  a  Beckman  and  Whitley  Model  #770  streak  camera  writing  at  J6mm/  s, 
and  a  shock  velocity  in  the  sample  calculated.  The  impedance  matching 
technique  was  used  to  calculate  the  corresponding  pressure  and  particle 
velocity  in  the  sample.  This  is  illustrated  in  Figure  2. 

Since  It  was  anticipated  that  the  M-30  propellant  might  react  during  these 
measurements,  an  additional  pressure  monitor  (a  Plexiglas  pellet)  was  placed 
on  the  propellant  free  surface.  The  shock  velocity  in  the  Plexiglas  pel  e 
was  measured  and  a  pressure  in  the  propellant  was  calculated  by  using  the 
interface  equation.  This  pressure  determination  in  conjunction  with  that 
obtained  using  the  impedance  matching  technique  was  used  to  insure^that  on  y 

unreacted  data  points  would  be  accepted  to  calculate  the  Hugcniot. 

If  the  propellant  were  to  react  It  would  produce  a  higher  pressure  in  the 
Plexiglas  overlay  than  would  unreacted  M-30.’  When  the  interface  equation  s 
used,  the  computed  particle  velocity  in  M-30  would  be  high  due  to  reaction  as 
illustrated  in  Figure  3A.  On  the  other  hand,  the  particle  velocity  calculated 
by  the  impedance  matching  technique  will  be  low  if  the  M-30  reacts  sufficien 
jy  to  increase  the  shock  velocity  In  the  sample. 

This  can  be  seen  by  referring  to  Figure  3B  where  the  particle  velocity  is 
determined  by  the  intersection  of  the  buffer  plate  release  adiabat  with  a  line 
having  a  slope  of  PQUg  where  Pq  is  the  initial  density  of  M-30  and  Ug  is  the 

shock  velocity  through  the  sample.  Thus  a  deviation  between  the  particle 


^Hice,  McQueen,'  Walsh,  "Compression  of  Solids  by  Strong  Shock  Waves,"  Solid 
State  Physics,  Volume  6,  1958.  , 

2Boyle,  Smothers,  Ervin,  "The  Shock  Hugoniot  of  Unreacted  Explosives,"  Fifth 
Symposium  on  Detonation,  August  18  -  21,  1970. 
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Figure  2.  Experimental  Arrangement  Used  to  Determine  the  Shock 
Velocity  in  the  Sample.  Also,  shown  is  a  graphical  illustration 
of  the  impedence-matching  technique  which  was  used  to  calculate 
the  pressure  and  particle  velocity  in  the  sample. 
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Figure  3A.  Graphical  Interpretation  of  the  Interface  Equation  for 
Reacting  and  Non-Reacting  M-30. 


Figure  3B.  Graphical  Interpretation  of  the  Impedence-Matching  Technique 
for  Reacting  and  Non-Reacting  M-30. 
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velocity  calculated  at  the  buffer-M-30  interface  and  that  calculated  at  the 
Plexiglas  -  M-30  in' crface  indicates  reaction.  Additional  information  is 
given  in  reference  2. 


IV.  RESULTS 

The  measured  data  points  for  M-30  propellant  are  shown  in  Table  1.  The 

numbers  in  parentheses  are  data  points  determined  using  a  Plexiglas  overlay 

and  the  interfare  equation  to  determine  particle  velocity  in  the  propellant. 

This  was  done  to  see  if  propellant  reaction  was  affecting  the  data  during  the 

t-.me  of  measurement  (less  than  one  microsecond).  The  first  four  lines  of  data 

showed  little  or  no  reaction  and  these  eight  points  were  used  to  determine  the 

unreacted  Hugoniot  for  solid  M-30  propellant.  A  least  squares  fit  of  the  . 

form,  U  =  b  *  su^  gave 
3  p 

U  =  2.661  X  105  +  1.655  u  cm/s. 
s  P 

The  measured  data  points  for  the  M-30  inert  simulator  are'  given  in  Table  2.  A 
least  squares  fit  to  these  data  points  gave  the  following  Hugoniot  for  the 
solid  inert  simulator. 

U  =  2.251  X  105  +  1.509  u  cm/s 
s  P 


Table  1.  M-30  Data. 


OS 

(gm/c  n^) 

Us 

(cm/s) 

UP 

(cm/sec) 

P 

(Kbars)* 

Comments 

1 .63 

3.515. X  105 

.508  1.542)  X  105 

29.6 

(31.6) 

1.66 

4. '103  X  105 

1.032  (1.081)  X  105 

75.4 

(79.0) 

1.66 

4.426  X  105 

1.031  (1.017)  X  105 

75.8 

(74.7) 

1.66 

4.391  X  105 

1.032  (1.118)  X  105 

75.3 

(81.5) 

1.66 

5.248  X  105 

1.464  (2.415)  X  105 

127.5 

(210.4) 

Reaction 

1.66 

5.226  X  105 

1.465  (2.258)  X  105 

127.1 

(195.9) 

Reaction 

1.66 

5.476  X  105 

Q 

1.445  (2.622)  X  ID5 

2 

131.3 

(226.1) 

Reaction 

*  1  Kilobar  =  10'  dynes/cm 
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Table  2.  M-30  Inert  Simulator  Data. 


3 

(go/cm  ) 

(cm/s) 

(cm/sec) 

(Kbars)* 

1.55 

3.050 

X  105 

.522 

X 

105 

24.7 

1.55 

3.051 

X  105 

.522 

X 

105 

24.7 

1.55 

3.863 

X  105 

1.064 

X 

105 

63.7 

1.55 

3.804 

X  105 

1.066 

X 

105 

62.9 

1.55 

4.611 

X  105 

1.548 

X 

105 

110.6 

1.55 

4.592 

X  105 

1.550 

X 

105 

110.3 

»  1  Kilobar  =  10^  dynes/cra2 

3  4  5  6 

The  porous  Hugoniots  were  calculated  from  the  solid  Hugoniots  using 

the  Mie-Gruneisen  equation  with  a  Gruneisen  gamma, 

calculated  from  thermodynamic  properties  of  the 


material  at  standard  temperature  and  pressure.,  At  standard  conditions,  the 

Gruneisen  gamma,  G  =  aC^/C^  where  a  is  the  volume  coefficient  of  thermal 

expansion,  Cq  is  the  bulk  sound  velocity  and  C  is  the  specific  heat. 

These  quantities  are  shown  in  Table  III.  The  volume  coefficient  of  thermal 
expansion  for  M-30  was  assumed  to  be 
-4 

3  X  10  cc  '  .  This  agrees  well  with  the  values  cited  in  the  literature 

co-°c  and  is  probably  a  good  approximation.  We  assumed  the  same 


value  for  the  inert  simulator.  The  bulk  sound  velocity  was  obtained  from  the 
unreacted  Hugoniot  of  the  solid  at  zero  particle  velocity.  The  specific  heat 
of  M-30  was  taken  from  the  literature.  The  specific  heat  of  the  inert  M-30 
simulator  was  measured  by  Leon  Decker,  BRL,  using  differential  scanning 
calorimetry. 


3 

Herrmann  W.,  "Constitutive  Equation  for  the  Dynamic  Compaction  of  Ductile 
Porous  Materials,"  JAP,  Vol.  40,  No.  6,  May  1969* 

4 

Hoffman,  Andrews,  Maxwell,  "Computed  Shock  Response  of  Porous  Aluminum,"  JAP, 
Vol.  39.  No.  10,  Sept.  1968. 

5 

Erkman,  Edwards,  "Computed  and  Experimental  Hugoniots  for  Unreacted  Porous 
High  Explosives,"  Sixth  Symposium  on  Detonation,  Aug.  24  -  27,  1976. 

^Zeldovich,  Raizer.  "Physics  of  Shock  Waves  and  High  Temperature 
Hydrodynamic  Phenomena,"  Vol.  II,  Academic  Press,  New  York,  1967. 


Table  3«  Grunei3en  Gamma. 


Material 

a 

Co 

CP 

G 

cc 

, (cm/sec) 

ergs 

0 

0 

cc-  c 

gm-  e 

M30 

.4 

3  X  10 

2.661  X  105 

1.502  X  107 

1.412 

Inert 

-4 

3  X  10 

2.251  X  105 

1.247  X, 107 

1.219 

The  porous  Hugoni  cs  were  calculated  using  the  following  simplifying 
assumptions: 

1.  G  remains  constant  at  its  zero  pressure  value  over  the 
range  of  interest. 

2.,  The  zero  pressure  specific  internal  energy  of  the  solid 
and  porous  material  are  equal. 

3.  The  materials  possess  zero  strength. 

4.  Mechanical  equilibrium  has  been  established. 

5..  Thermal  equilibrium  has  been  established.  Refer  to  the. 

Discussion  Section  for  more  on  assumptions  4  and  5. 

The  Mie-Gruneisen  equation  of  state  relates  the  thermal  component  of  pressure 
to  the  thermal  component  of  specific  internal  energy  at  a  given  volume.  This 
can  be  expressed  by  the  following  equation: 

<pf  -p3)  .0  (Ef-E.)  O) 

where  the  subscript  f  refers  to  the  porous  material  and  the  subscript  s  refers 
to  the  solid  material.  The  validity  of  this  equation  is  treated  in  references 
3,  4,  5  and  6. 


Using  the  Rankine-Hugoniot  energy  equation  for  the  porous  and  the  solid 
materials  the  following  equations  can  be  written; 


Ef  -  Eof  = 

1/2 

Pf 

U-i 

o 

> 

-  V) 

Es  -  Eos  * 

1/2 

Ps 

(Vos 

-  V) 

(2) 

(3) 


Refer  to  Figure  4  for  a  geometrical  interpretation  of  equations  (1),  (2),  and 
(3;.  Since  we  assumed  E  ~  =  E  equations  (1),  (2),  and  (3)  can  be  combined 
to  give:  01  03 


(4) 


Since  the  Hugoniot  of  the  solid  material  has  been  determined 
experimentally  (Ug  =  b  +  su  ),  we  can  use  this  relationship  and  the  equations 

for  conservation  of  mass  and  momentum  across  the  shock  wave  to  derive  an 
expression  for  P  ,  the  shock  pressure  in  the  solid  material, 


•  *  1  -  "  1  »  J  /  V  y  ’*  i  ay  ~ 


( 


SOLID  HUGONIOT 


POROUS  HUGONIO' 


Figure  4.  Geometrical  Interpretation  of  Equations  (1),  (2),  and  v 3 ) . 

(pt  -  Va  ■  \  pf  <vof  ■  v)  •  j  p8  <v°=  •  v)  ■ 


Equation  (4)  can  then  be  rewritten, 


b2  -  v) 

Ks  ^os-W 


(6) 


P^.,  the  pressure  in  the  porous  material,  is  then  calculated  by 
substituting  values  for  V,  the  spacific  volume  of  the  shocked  material, 
into  Equation  (6).  i'f  will  be  in  kilobars  for  V  in  emVgm  and  b  in  cm/s. 
Conservation  of  mass  across  the  shock  front  requires 


Usf  .  Usf  -  %f 

Conservation  of  momentum  requires 


(7) 


The  porous  Hugoniot  values  calculated  using  equations  6,  9,  and  10  are  shown 
in  Tables  4  and  5. 


Table  4.  Porous  Inert  Hugoniot  Calculations* 


G 

V 

V 

V 

U3f 

upf 

(cm^ /gm) 

3 

(cm  /gm) 

(Kbar) 

(cm/s) 

(cm/sec) 

1.216 

1.250 

.63 

4.9 

1.109  X 

105 

.550  X  105 

1.216 

1.250 

.62 

8.8 

X 

r— 

CO 

•=JV 

• 

105 

.747  X  105 

1.216 

1.250 

.61 

13.6 

1.821  X 

105 

.932  X  105 

1.216 

1.250 

.60 

19.3 

2.152  X 

105 

1.119  X  105 

1.216 

1.250 

.59 

26.2 

2.490  X 

105 

1.315  X  105 

1.216 

1.250 

.55 

34.7 

2.844  X 

105 

1.525  X  105 

1.216 

1.250 

.57 

45.3 

3.226  X 

105 

1.755  X  105 

1.216 

1*250 

.56 

58.8 

3.648  X 

105 

2.014  X  105 

1.216 

1.250 

.55 

76.? 

4.125  X 

105 

2.310  X  105 

1.216 

1.250 

.54 

99.6 

4.681  X 

105 

2.659  X  105 

1.216 

1.250 

.53 

131.9 

5.351  X 

105 

3.082  X  105 

1.216 

1.250 

•  52 

179.2 

6.193  X 

105 

3.617  X  105 

1.216 

1.250 

.51 

253.6 

7.137  X 

105 

4.332  X  105 

1.216 

1.250 

.50 

385.5 

8.961  X 

105 

5.337  X  105 

A  least  squares  fit  of  the  shock  and  particle  velocities  shown  above  gives 
the  porous  inert  Hugoniot. 


Table  5.  Porous  M-30  Hugonlot  Calculations, 


G 

Vof  , 

V 

Pf 

Usf 

V 

(cm  /gm) 

3 

(cm  /gm) 

(Kbar) 

(cm/s) 

(cm/s) 

1.412 

1.149 

.59 

7.7 

1.351  X  105 

.658  X  105 

1.412 

1.149 

.58 

15.7 

1.911  X  105 

.947  X  105 

1.412 

1.149 

.57 

25.9 

2.431  X  105 

1.225  X  105 

1.412 

1:149 

.56 

39.1 

2.959  X  105 

1.517  X  105 

1.412 

1.149 

.55 

56.5 

3.528  X  1,05 

1.839  X  105 

1.412 

1.149 

.54 

80.2 

4.169  X  10S 

2.210  X  105 

1.412 

1.149 

.53 

113.8 

4.926  X  105 

2.654  X  105 

1.412 

1.149 

.52 

164.1 

5.869  X  105 

3.214  X  105 

1.412 

1.149 

.51 

245.9 

7.129  X  105 

3.966  X  105 

1.412 

1.149 

.50 

398.7 

9.006  X  105 

5.088  X  105 

A  least  squares  fit  of  the 
porous  M-30  Hugoniot. 

shock  and  particle 

velocities  shown 

above  gives  the 

The  Hugoniot 
by: 

of  porous 

M-30  (bulk  density 

■j 

=  .870  gm/cm  )  is 

represented 

°.f  ■  -3’8  x 

105  +  1.720  u  .  cm/s. 
Pf 

The  Hugoniot  of 
by: 

the  porous 

inert  (bulk  density 

3 

=  *800  gm/cm  )  is  represented 

°af  *  *337  X 

5 

10  +  1.617  u  .  cm/sec, 

Pf 

These  equations 

are  least  i 

squares  fit3 

to  the 

calculated  porous 

Hugoniot 

values.  There  is  no  physical  significance  to  the  intercept  value  at  zero 


particle  velocity.  A  power  function  would  provide  a  better  fit  to  the  low 
velocity  points.  The  Hugoniots  in  the  pressure  -  particle  velocity  plane  are 
shown  in  Figure  5. 


(kbar) 


V.  DISCUSSION 


The  porous  Hugoniot  data  allow  one  to  calculate  the  pressure  in  the 
propellant  bed  due  to  shaped  charge  jet  impact.  This  information  is  useful  in 
performing  hydrodynamic  calculations  for  a  non-reacting  propellant  bed.  As 
shown  in  Table  1  the  solid  propellant  undergoes  violent  reaction  at  an 
incident  pressure  of  approximately  127  kbars,  increasing  to  around  200  kbars 
within  the  measurement  interval  of  one  microsecond.  However,  it  is 
anticipated  that  reaction  could  occur  at  lower  pressures  of  longer  duration, 
conditions  which  exist  in  the  pressure  field  generated  by  jet  impact.  The 
solid  Hugoniot  reaction  pressure  may  be  a  good  way  to,  a  priori,  compare  the  , 
response  of  various  propellants  upon  shaped  charge  jet  impact  since  impact 
pressures  can  be  determined  from  known  jet  velocity  and  Hugoniots  for  the  Jet 
and  propellant  bed. 

The  question  of  the  approach  to  mechanical  equilibrium  of  porous  samples 
is  treated  in  reference  4.  This  reference  used  a  one-dimensional  Lagrangian 
computer  program  to  simulate  the  response  of  powdered  aluminum  to  a  shock. 

The  powdered  aluminum  was  mocked  up  by  a  series  of  flat  plates  separated  by 
gaps  of  width  equal  to  the  plate  thickness.  The  calculated  results  were 
compared  to  experimental  results  on  powdered  aluminum  and  found  to  be  in  good 
agreement.  The  calculations  showed  an  approach  to  equilibrium  behind  the 
shock  for  porous  aluminum;  the  equilibrium  states  were  consistent  with  the 
Rankine-Hugoniot  jump  conditions  applied  to  the  porou  aluminum.  The 
computational  results  indicated  that  tne  shock  had  to  travel  about  ten  times 
the  plate  thickness  before  mechanical  equilibrium  was  established  behind  the 
shock  wave. 

The  three-dimensional  nature  of  a  propellant  bed  would  tend  to  promote  a 
more  rapid  equilibration  by  the  increased  number  of  shock  interactions  per 
unit  particle  per  unit  time.  On  the  other  hand,  "hot  spots"  would  be  produced 
in  regions  of  convergent  mass  flow  and  because  of  the  large  grain  size,  these 
"hot  spots"  could  persist  over  relatively  long  times.  Therefore,  the 
assumption  of  thermal  equilibrium  may  not  be  strictly  valid.  Despite  this 
limitation.  It  is  believed  that  the  porous  Hugoniots  calculated  here  are 
useful  approximations  to  a  real  propellant  bed. 
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LIST  OF  SYMBOLS 


SYMBOL 

a  volume  coefficient  of  thermal  expansion 

b  coefficient  i.n  the  Hugoniot  relationship,  U  =  b  +  su 

s  p 

3  coefficient  in  the  Hugoniot  relationship,  U  =  b  +  su 

a  p 

C  bulk  sound  velocity 

o. 

Cp  specific  heat  at  constant  pressure 

E  specific  internal  energy 

E  specific  internal  energy  of  compressed  porous  material 

E  „  initial  specific  internal  energy  of  porous  material 
of 

£g  specific  internal  energy  of  compressed  solid  material 
Eq3  initial  specific  internal  energy  of  solid  material 
G  Gruneisen  constant 

P  pressure 

Pp  pressure  in  compressed  porous  material 

PQf  initial  pressure  in  porous  material 

Pg  pressure  In  compressed  solid  material 

P  initial  pressure  in  solid  material 

os 

U  3hock  velocity 

U  shook  velocity  in  porous  material 

Up  particle  velocity 

u  particle  velocity  in  porous  material 

V  specific  volume  of  compressed  material 

VQf>  initial  specific  volume  of  porous  material 

V  initial  specific  volume  of  solid  material 

os 

P  density  of  compressed  material 

p  initial  density  of  solid  material 

os 

p  .  initial  density  of  porous  material 
of 
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